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We find that collective flow model which can successfully analyze charged particle distributions 
at AGS and lower SPS (-Ejvat l ess than 20GeV in the lab frame), but fails to analyze that of at 
RHIC. The tails of distribution of charged particle at RHIC has a jump from the collective flow 
model calculation as the energy increases. Thermalization Component Model is presented based on 
collective flow to study the multiplicity distributions at RHIC in this paper. It is realized that the 
region of phase space of collective flow can reflect that of thermalization region. By comparing the 
contributions of particle productions from thermalization region at different energies and different 
centralities, we can deepen our study on the feature of collective movement at RHIC. 
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I. INTRODUCTION 

One of the central question at RHIC is the extent 
to which the quanta produced in the collision interact 
and thermalize [l|, [|J ■ Nuclear collisions generate enor- 
mous multiplicity and transverse energy, but in what ex- 
tent does the collision generate matter in local equilib- 
rium which can be characterized by the thermodynamic 
parameters temperature, pressure, and energy density? 
Only if thermalization has been established can more de- 
tailed questions be asked about the equation of state of 
the matter. 

Recently it is realized that the study of collective flow 
is one of the important tools to study multi-hadron pro- 
duction of relativistic heavy- ion collisions @, H, ■ This 
is because the longitudinal and transverse flow includes 
rich physics, and collective flow relates closely to early 
evolution and nuclear stopping. Collective flow is often 
utilized to express the thermalization degree of relativis- 
tic heavy-ion collisions system. Detailed studies of the 
observed final state flow pattern will deep our under- 
standing of dynamic mechanism of relativistic heavy-ion 
collisions. Collective Flow Model HH (CFM) was devel- 
oped basing on the pure thermal model. It achieves suc- 
cess in the discussion of charged particle distribution at 
AGS and lower SPS energy (20GeV) and become an indi- 
cator for the existence of collective flow at AGS 0,111]. But 
a detailed analysis of the experimental data at SPS and 
RHIC energy with CFM has shown that as the increase 
of collision energy, the two tails of the charged hadron 
distributions have a symmetric jump away from the cal- 
culation of CFM. These phenomena will make us to re- 
consider collective flow theory at higher collision energy. 

As shown in Fig.l, CFM fails to analyze the charged 
hadron distribution at RHIC energy regions. As the in- 
crease of the collision energy, the tail of the distribution 
of experimental data jump from the calculation of CFM. 



The naive reason seems to be that experimental data on 
hadron yields are now available over a broad collision 
energy range as the increase of collision energy. As the 
increase of phase space of particle distribution, thermal- 
ization at whole phase space of particle production be- 
comes more difficult. Detailed analysis of thermalization 
relation with centrality and energies at RHIC is needed. 

How to simulate the data of charged hadron distribu- 
tion at higher SPS and RHIC energy regions and what 
these results tell us are the main topic of the paper. PHO- 
BOS has used three Gaussian distributions to simulate 
the distribution of charged hadron successfully [8l4lCt Il3j . 
Georg Wolschin [TJ] et al also discussed the charged 
hadron distribution by using three component distribu- 
tion functions to construct FokkerPlank equation. Their 
models both assumed three random Gaussian distribu- 
tion emitting sources. The Three Gaussian sources rep- 
resent target, projectile and central source in physics, 
respectively. 

The main goal of this paper is to give a study of the 
thermalization features of multi-particle production in 
heavy ion collision at high energy in the framework of 
the collective flow theory. We restrict ourselves here to 
the basic features and essential results of the CFM ap- 
proach. A complete survey of the assumptions and re- 
sults, as well as of the relevant references, is available in 
Ref.di, [H-[I1. 

The paper is organized as follows. The analysis details 
based on the Thermalization Component Model(TCM) are 
described in Sec. 2. The comparisons of TCM calculations 
with experimental data and the related theoretical anal- 
ysis with TCMare given in Sec. 3. A summary is given in 
Sec. 4. 
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FIG. 1: (a) 7r meson rapidity distribution of £^=30 GeV 
at SPS [3; (b) charged hadron pseudo-rapidity distribution 
of 75^7=200 GeV 



II. THERMALIZATION COMPONENT MODEL 

The hot and dense matter produced in relativistic 
heavy ion collisions may evolve through the following sce- 
nario: pre-equilibrium, thermal (or chemical) equilibrium 
of partons, possible formation QGP or a QGP hadron gas 
mixed state, a gas of hot interacting hadrons, and finally, 
a freeze-out state when the produced hadrons no longer 
strongly interact with each other. Since the produced 
hadrons carry information about the collision dynamics 
and the entire space-time evolution of the system from 
the initial to the final stage of collisions, a precise anal- 
ysis of the multiplicity distributions of charged hadrons 
is essential for the understanding of the dynamics and 
properties of the created matter. 

A detailed analysis of the experimental data at 
SPS and RHIC energy with CFM has shown that as 
the increase of collision energythe two tails of ir or the 
charged hadron distributions show a (symmetric) dis- 
crepancy between the data and the calculation. These 
phenomena will make us to reconsider Collective flow the- 
ory at higher collision energy. Detailed analysis of the re- 
lation with thermalization with centralities and energies 
at RHIC is needed. Let us first sketch our overall picture 
and detail our arguments subsequently. The model we 
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FIG. 2: The pseudo-rapidity distributions at ^/smn =19.6, 
62.4, 130, 200 GeV for Au+Au Collisions. Experimental data 
are given by triangle [jj-flol].The solid lines are the results 
given by TCM, which is the summation of the three compo- 
nent contributions. 



considered contains three distinct assumptions some of 
which arc rather different from those usually contained 
in other flow models. 

(i) The size of phase space of the particle distribution 
increases with the increase of collision energy. It seems 
more difficult to realize thermalization at the whole phase 
space of particle production at SPS and RHIC data. It 
is assumed that the Gaussian distributions were fit to 
the distributions of the produced charged hadrons at the 
two fragmentation regions, and thermalization prefers to 
occur at the central rapidity region at SPS and RHIC. 

(ii) The collective flow of central rapidity region carries 
information of the early time of heavy-ion collision. The 
system expands not only in the longitudinal direction, 
but also in the transverse direction. The two dimensional 
collective flow is used to study the thermalization process 
at RHIC. 

(iii) The phase space is compartmentalized as the ther- 
malization region and non-thermalization regions. The 
non-thermalization regions locate at the two fragmenta- 
tion regions. The total multiplicity distributions are the 
summation of the contributions from the target fragmen- 
tation region, projectile fragmentation region and central 
region, respectively. 
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Here i = 1,2,3 denotes target, projectile and central 
region, respectively. Ni and Fi are the particle numbers 
and the normalization functions of target, projectile and 
central regions, respectively. 
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As assumed before, the distributions of target and pro- 
jectile fragmentation regions are given with Gaussian dis- 
tributions: 
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Here a is the distribution width of Gaussian, y\ , y 2 are 
the locations of central of target and projectile emitting 
source. 

F3 is the distribution of two dimensional flow, which 
is given by 0,0] 
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Here m l t ° and mj" are the experimental limits in which 
the spectrum is measured. The freeze-out radius Rf and 
the longitudinal extend of the fireball is fixed via the 
finite interval (-770, 770), lo'is modified Bessel function. 

For the two dimensional flow theories, we should say 
a few words. The geometry of the freeze-out of two di- 
mensional flow hyper-surface crj fixed as follows: in the 
time direction we take a surface of constant proper time 
. In r)i direction, the freeze-out volume extends only to 
a maximum space-time rapidity 770 , which is required 
by the finite available total energy and breaks longitudi- 
nal boost- invariance proposed by Bjorken [l6| . In the 
transverse direction the boundary is given by Rf, which 
describes a cylindrical fireball in the 77 — r space. The 
detailed discussion was shown in Ref. 0, [H] ■ 



III. COMPARING WITH EXPERIMENTAL 
DATA 

It is found that CFM describe experimental data of 
charged particle distribution very well when we discuss 
Au-Au center collisions at AGS energy region. The con- 
tribution of fragmentation regions can be ignored, so ex- 
pressions (I) can be predigested: 
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The results from CFM are consistent with experimen- 
tal data in Au-Au collisions at AGS energy region, such 
as E l £ b N =2, 4, 6, 8, 11.6 GeV in the lab frame. This in- 
dicates that when at lower AGS energy region CFM can 
describe the charged particle distribution well, then our 



thermalization component model revert to collective flow 
model. The reason seems that phase space is small and 
the nucleus stopping power is very strong at AGS energy 
region, so particles can be almost completely thermalized 
in whole phase space. The same situation is true for that 
of SPS energy region below 20 GeV. 

But with the collision energies increase (E l fi b N above 30 
GeV), the experimental points have a symmetric jump 
away from the calculation of CFM at two tails (as shown 
in Fig.l). This phenomenon can be explained by the nu- 
clei's penetrability. The higher the collision energies, the 
more transparent the nuclei, and the larger the extension 
of the phase space of the produced particle. Collective 
flow is formed at the central rapidity region after thermal- 
ization. The distributions of non-thermalization charged 
hadrons are presented by Gaussian. The thermalization 
area becomes one part of the whole phase space. 

Since June 2000, the Relativistic Heavy-Ion Collider 
(RHIC) has opened a new energy region for the study of 
multi-hardon production. We have analyzed the exper- 
imental data of charged particle distribution in Au-Au 
center collisions in the RHIC energy region from 19.6 to 
200 GeV of y/I 

NN 

We can calculate the rapidity distribution of charged 
particles whit expressions (1). It is known that we trans- 
fer rapidity distribution to pseudo-rapidity distribution 
just by multiplying a factor [2~i| : 
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FIG. 3: The dependence of the percentage of the charged 
hadron production from the thermalization regions on the col- 
lision energies. 
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FIG. 4: The charged hadron pseudo-rapidity distribution 
at different centrality at y / iiv"iv=62.4GeV and 200GeV for 
Au+Au and Cu+Cu Collisions, respectively. Solid lines are 
the results from TCM, Experimental data are given by PHO- 
BOS l-fil 



We fit the experimental data of RHIC energy region by 
HCM model by x 2 /dof . The comparison of the measured 
and calculated distributions for the best fit ( x 2 /dof min- 
imization) is presented in Fig. 2. The TCM calculations 
are accordant with experimental data shown in Fig. 2. 
The percentages of the charged hadron productions from 
the thermalization regions at AGS, SPS and RHIC energy 
region are presented in Fig. 3 by TCM. It is found that 
most of the produced particles at AGS come from the 
thermalization region, and the percentage of produced 
particles from the thermalization region decreases as the 
energy increase. The reduction trend becomes weaker 
and seems to reach saturation as ^Jsnn reaches 62.4GeV 
at RHIC energy region. The detailed fit parameters of 
our TCM with experimental data are shown in Table 1. 



TABLE I: The fit results of TCM with the experimental data 
at SPS and RHIC energy regions 
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collisions at snn =62.4 and 200 GeV. It is found that 
the calculation results from TCM are consistent with that 
of the experimental data. The results are presented by 
Fig. 4 and Table 2. The experimental data are taken from 
Ref. 
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PHOBOS Collaboration Working at RHIChas presented 
many experimental data [8l-fl0| of different energy and 
different centrality including Au-Au collisions and Cu-Cu 
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FIG. 5: The dependence of the percentage from the thermal- 
ization region on different centralities for ^snn =62.4, 200 
GeV 

It is shown from Fig. 5 that the percentage ratios of the 
particle production from the thermalization regions in- 
crease with the increase of the centralities at RHIC. From 
Fig. 5 (a), It is found that the contribution ratios from the 
thermalization region is appreciably larger for the smaller 
collision system (Cu ±Cu) than that of larger collision 
system (Au ±Au) at ^/snm— 62.4GeV. But from Fig. 5 
(b), we find that the percentage ratios of particle produc- 
tion from thermalization regions is almost independent of 
the size of collision systems at ^/snn— 200 GeV . 

In our TCM, the free parameters are the limitation of 
collective flow r/o and the emission sources' positions in 
fragmentation area 2/1,2- We have y\ = —yi in the case 
of symmetry collisions. The values of transverse flow 
and temperature of collective flow refer to Ref. [1, [H UH, 
[20j . The values of rii(i = 1, 2, 3) are numbers of particles 
from the fragmentation and the thermalization regions, 
respectively. 

A linear relationship is obtained between 770 and 
In^/sNN by detailed study. The linear equations are 
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given by fitting four data at SPS 
as follows: 



RHIC energy regions 



r/o = OAOIut/s^ + 0.71 



(7) 



here 770 is the extension of collective flow. From Eq.7, 
we can predict the extension of the thermalization region 
at LHC with the collision energy increase. 
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FIG. 6: The relation between the limitation of thermalization 
region with Iu^/smn ■ 

Here, we should mention that quite a few theoreti- 
cal models can give equally good representation of the 
data of particle productions at AGS,SPSand RHIC, such as 
these thermal models (22h25| based on the assumption of 
global thermal and chemical equilibrium, and hydrody- 
namic models [26l - [3l| based only on the assumption of lo- 
cal thermal equilibrium, to transport models [32h37| that 
treat nonequilibrium dynamics explicitly. The thermal 
models have been very successful in accounting for the 
yield of various particles and their ratios, while the hydro- 
dynamic models are particularly useful for understand- 
ing the collective behavior of low transverse momentum 
particles such as the elliptic flow. Since transport mod- 
els treat chemical and thermal freeze-out dynamically, 
they are also natural and powerful tools for studying the 
Hanbury-Brown-Twiss interferometry of hadrons. 

For hard processes that involve large momentum trans- 
fer, approaches based on the perturbative quantum chro- 
modynamics (pQCD) using parton distribution functions 
in the colliding nuclei have been used [38]. Also, the clas- 
sical Yang-Mills theory has been developed to address 
the evolution of parton distribution functions in nuclei at 
ultra-relativistic energies and used to study the hadron 



rapidity distribution and its centrality dependence at 
RHIC|39 41] . These problems have also been studied in 
the pQCD-based final-state saturation model [42Tj4^ ]. A 
multiphase transport (AMPT) model that includes both 
initial partonic and final hadronic interactions and the 
transition between these two phases of matter (45M48j | was 
constructed to describe nuclear collisions ranging from p 
- A to A - A systems at center-of-mass energies from 
about ^/sWn = 5 to 5500 GeV at LHC. 



IV. SUMMARY AND CONCLUSIONS 

Hadron multiplicities and their distributions are ob- 
servables which can provide information on the nature, 
composition, and size of the medium from which they are 
originating. Of particular interest is the extent to which 
the measured particle yields are showing thermalization. 
The feature of thermalization of high energy heavy ion 
collisions at RHIC has been analyzed in this paper. 

CFM fails to analyze the charged particle distributions 
when the collision energies increase to above 30 GeV. The 
tail of distribution of the charged particle at RHIC has a 
jump from the CFM calculation with the energy increase. 
The naive reason seems to be that the experimental data 
on hadron yields are now available over a broad colli- 
sion energy range with the increase of collision energy. 
It seems more difficult for thermalization at the whole 
phase space of particle production with the increase of 
the phase space of particle distribution. 

On the other hand, the phenomena may suggest that 
something else happens, including interaction mech- 
anism, such as the onset of de-confinement in the 
early stage of the reaction with the collision energy 
(-Etvw) above 30 GeV at the lab frame, which has been 



mentioned in Ref.[49(. In Ref.[49|, central Pb - Pb colli- 
sions were studied in the SPS energy range. At around 
-^nn ~30 GeV the ratio of strangeness to pion production 
shows a sharp maximum, the rate of increase of the pro- 
duced pion multiplicity per wounded nucleon increases 
and the effective temperature of pions and kaons levels 
to a constant value. These features are not reproduced 
by present hadronic models, however there is a natural 
explanation in a reaction scenario with the onset of de- 
confinement in the early stage of the reaction at SPS en- 
ergy. 

Collective flow in heavy-ion collisions is an unavoid- 
able consequence of thermalization. The extension of the 
phase space of collective flow can reflect that of thermal- 
ization region. It is found that the TCM can fit the 
experimental data well for the particle production at the 
whole AGS, SPS and RHIC energy regions. The per- 
centage ratios of contributions of the particle production 
from the thermalization region are the largest at AGS, 
and decrease as collision energies increase at SPS and 
RHIC, but seem to reach saturation when y / sjv r /v=62.4- 
200 GeV at RHIC. It is also found that the extension of 
the flow shows a linear dependence on Iu^/snn- From 
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that, we can predict the thermalization extension at fu- 
ture LHC experimental data. 

It is shown from our study that the percentage ratios of 
particle production from thermalization regions increase 
with the increase of the centralities at RHIC. The con- 
tribution ratios from thermalization region are apprecia- 
bly larger for the smaller collision system (Cu + Cu) at 
y / s]v]v=62.4GeV, but independent of the collision system 
at y/JWN=200 GeV . 
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